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1 ABSTRACT 

t 

: » 

i 

The  concept  of  a compliant  wall  wind  tujuiol , whoso 
wall  deforms  in  response  to  aerodynamic  loading,  is  studied 
theoretically.  It  is  argued  that  such  walls  have  advan- 
tages for  the  design  of  wind  tunnels  with  minimum  wall 
interference. 
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♦ 
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LIST  OF  SYMBOLS 


- see  equation  (17) 

~ see  equation  (25) 

- wind  tunnel  height 

- airfoil  chord 

- see  equation  (23) 


- torsional  spring  stiffness 

- aerodynamic  moment  about  hinge 

- free  stream  Mach  number 

- see  equation  (17) 

- aerodynamic  pressure 

- Laplace  transform  variable 

- see  equation  (17) 

- see  equation  (25) 

- free  stream  velocity 

- downwash 

- streamwise  coordinate 

- transverse  cordinate 

- angle  of  attaclc 
r (M^  - 1)^^^ 

= Poo  (eb)V6K 

- wall  deflection 

- velocity  potential 
~ free  steam  density 


~ compliant  wall  rotation 


I 


'•  ■--•V ■■  .r 


0^  - compliant  wall  rotation  corresponding  to  zero 


spring  twist  or  X = 0 


= 0-0. 


Superscripts 


* - Laplaco  transform 

TE  - trailing  edge 

ZI  - zero  interference 

% - renormalization  by  multiplying  or  dividing  by 


Subscript 


- divergence 


I' 


INTRODUCTION 


'I  The  self-correcting  wind  tunnel  hus  attracted  con- 

siderable interest  in  recent  years  as  a concept  for 

I 

M achieving  minimum  wall  interference  in  a wind  tunnel. 

I Sears,  et.  al.^,  have  presented  an  authoritative  account 

i which  is  reconunended  to  the  reader.  Among  their  important 

I ! 

[■)  conclusions  is  that  a perforated  wind  tunnel  wall  (with  con- 

i 

trol  of  the  pressures  in  subdivided  plenum  chambers)  is 

superior  to  control  of  impermeable,  flexible  walls  such  as 

I 2 I 

suggested  by  Cheva liter  and  Goodyer' . In  Ref.  2 and  3 and 

( 

the  discussion  in  Ref.l  it  is  assumed  that  tlie  imperme- 
able, flexible  wall  is  deformed  into  an  apriori  determined 
shape  of  fixed  magnitude.  The  selection  of  wall  shape  func- 
tion and  its  magnitude  is  made  so  as  to  minimize  wind 
tunnel  wall  interference.  In  the  present  paper  this  con- 
cept is  broadened  to  allow  at  least  the  magnitude  of  the 
wall  shape  to  vary  with  flow  conditions,  e.g.  flow  dynamic 


pressure.  More  generally  the  shape  function  might  vary  as  well. 
As  will  be  shown  there  are  advantages  to  such  a compliant 
wall  concept. 

For  simplicity  we  consider  two-dimensional,  supersonic 
flow  in  the  context  of  small  perturbation  theory.  The 
compliant  wall  is  modeled  as  rigid  hinged  segments  whose 
compliance  is  determined  by  attached  springs.  As  the  reader 
will  appreciate  the  basic  concept  allows  for  a much  wider 
range  of  flow  regimes  and  compliant  wall  characteristics 
at  the  cost  of  some  increased  complexity  in  mathematics  and, 
ultimately,  experimental  implementation. 


ANALYSIS 


The  fluid  equation  of  motion  for  two-dimensional , 
small  perturbat  ioi^ , potential  flow  is 


p if  - if  =0 

^xx  zz 

where  P = (M^,  - 1)^'^^ 


The  boundary  condition  on  the  airfoil  is 


whore  w is  the  airfoil  downwash.  For  an  airfoil  at 
steady  anqlc  of  attack,  a,  w = - The  bovindviry  con- 

dition on  the  wind  tunnel  wall  is 


= 11  n 

X 


z=b/. 


where  b is  the  wind  tunnel  heiqht  and  n the  deflec- 
tion shape  of  the  wind  tunnel  wall. 

To  solve  (1),  subject  to  (2)  and  (S),  a Laplace 
Transform  with  respect  to  x is  used.*  Solvinvj  the  re- 
sulting ordinary  differential  equation  in  z and  using 
the  tranformed  boundary  conditions  gives  tor  the  l.a- 

place  transform  of  the  velocity  potential,  if*, 

T ,w*  s Pb/-  ,,  *1  * 

If*  = 2 (-  e 2 - U^,n*l  ^ w* 

2=0  ^ (4) 

e®  ^^'^2  - e‘®  ^^^^2 


♦following  Miles,  Ref. 4. 
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+ U 1 cot  h ^ 


z-b/. 


(^s 


+— 
{^S 


lo 


(5) 


Our  primary  interest  is  in  the  fluid  pressure,  p,  which 
is  determined  from  Bernoulli's  equation. 


P 


vX» 


u 


4' 


X 


(6) 


To  facilit.ate  the  invt'rsions  of  (-J)  and  (TO  and  sub- 
sequently determine  p from  (b) , the  followinq  identi- 
ties are  used. 


-st'b/^  V -nst'b 
= e 2 : e 


(V) 


n=0 


coth  st'b/^  = 1 + 2^;  e 
^ n=l 


-nst'b 


(8) 


Using  these  identities  and  standard  inversion 
formulae,  the  following  results  are  obtained. 


z=0 


^ f w (x)  + 2 i;  , , ... 

= - ^ ^ I-  w(x-n  Bb) 


+ U r,  (X-  eb  -(n-l)t'b)  1 ) 

X 2 


(9) 


p 


m^b) 


z=b/^ 


uo 


P u 

CW  CO 


p 


}: 

n=l 


w (x+ 


2 


P 


In  (x)+2  X n (x-nPb) 1 
^ n=l 


(10) 


VVo  shall  use  (10)  subsoqiiently ; first  wo  turn  our  attention 
to  (9).  The  first  term  in  (9)  is  the  value  of  p if  there 
were  no  wind  tunnel  wall  (interference);  the  second  term 
(first  sumnivition)  is  the  effect  of  a riejid  hori;^ontal  wall; 
the  third  term  (second  summation)  is  the  effect  of  the  wind 
tunnel  wall  deflection  sliape  with  respect  to  the  horizontal 
reference  line.  For  zero  wall  interference  we  choose  ti 
such  that  the  second  and  third  terms  cancel.* 

Thus  ' 

n^^(x)  = g (x-pb/^)  (11) 


One  way  of  satisfying  (11)  is  to  have  a rigid  but  moveable 

wall  (perhaps  composed  of  hinged  segments)  which  can  be 

controlled  perfectly  to  obtain  an  appropriate  n.  This  is  cssen- 

o 3 

tially  the  concept  of  Chovallier  and  Ooodyer  . However,  any 
practical  moveable  walls  will  have  some  compliance  or  flexi- 
bility; moreover  as  wo  shall  show  there  may  be  some  advantages 
to  constructing  a deliberately  compliant  wall. 

SINGLE  SEGMENT  WITH  HINGE  AT  LEADING  EDGE: 

For  definiteness,  we  consider  a zero  thickness  airfoil 

* Note  there  is  no  wall  interference  if  c^pb  since  t)ien  the  second 

and  third  terms  arc  identically  zero.  Physically  the  waves  reflected 
by  the  wall  do  not  strike  the  airfoil. 
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at  constant  angle  of  attack,  a,  and  a compliant  wall  con- 
sisting of  a single  rigid  segment  hinged  at  its  leading  edge. 
The  compliance  is  due  to  an  elastic  torsional  spring  as 
shown  in  Fig.  2.  It  is  assumed  the  leading  edge  of  the  wall 
segment  begins  at  x = in  practice  this  might  require  a 

sliding  leading  edge  to  accomodate  changing  p or  a multi - 
segment  wall,  each  segment  of  which  can  be  locked  out  or  in 
as  needed.  The  wall  segment  has  its  trailing  edge  at  x = 

c + P1V2_ 

The  wall  angle  for  zero  spring  stretching  is  0^;  the 
spring  twist  angle  is  AO.  Thus  the  total  wall  segment  angle, 
0,  is  given  by 

0 = 0^  + AO  (12) 

For  the  geometry  in  question, 

\ 

w = - Ujj,  n for  0 < X < c 

n = - 0 for  pb<  X < pb  + c (13) 


The  equation  of  equilibrium  for  wall  segment  is 


M + K (0  - Oq)  = 0 


(14) 


where  K is  the  torsional  spring  constant  and  M the  aerody- 
namic moment  about  the  hinge  point  (positive  in  the  direction 
of  rotating  the  segment  to  close) . M is  given  by 


M =/ 


c+ 


Pb 


Pb/. 


z=b/. 


(x-pb/2 )dx 


(15) 


Substitution  of  (13)  into  (10)  and  the  result  into  (15) 
gives 
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:-z  , 


M = - K Xlu  - 0 Tjjl 


(16) 


whore 


X = p U ^ (pb)^,  is  an  aeroclastic  coefficient. 


6 K 
N-1 

A„  5 ):  n (2n-l ) 

n=l  2 


+ N - (N-1) ] 


c - (N-1) 

Bb 


+ (N-1)]  , 


(N-1)  < p b 


(17) 


Substitution  of  (16)  into  (14)  and  solving  for  0 gives 


0 

a 


(18) 


N 


a 


1 + XT. 


H 


N-1 

T = K 
^ n=l 


is  an  aerodynamic  coefficient,  and 


(2n^l) 

2 


+ (2N-1)  (N-l)l 

• l|j^-  (N-1)  +(N-1)1  , 

2 

is  an  aerodynamic  coefficient,  when'  tlie  integer,  N,  is  determiiu'd  by 
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HINGE  AT  TRAILING  EDGE: 

A similar  calculation  for  the  wall  seyment  hinged  at 
its  trailing  edge  rather  than  its  leading  edge  gives 


a n.'l'E  J 

K1 


(19) 


1 - AT 


tE 


Note  that  aeroelastic  divergence  occurs,  0/  ► «>,  wlien 

CONDl'I’IONS  FOR  ZERO  WAJJ-  INTERFERENCE: 

From  (11)  and  (13),  for  zero  wall  interference  we  require 
(0/^)^'^  = 1 (20) 

Using  (20)  and  (18)  or  (19),  the  aeroelastic  coef f ic i I'n t at 
which  zero  wall  interference  occurs,  A , may  In'  comj'Uted  as 

for  hinge  at  the  leading  edge  (21) 


0, 


u 


or 


, ZI 


T - A 
N N 


{■  ■ 

TE  _ TE 
N N 


for  hinge  at  the  trailing  edge  (22) 


Since  A^'^  must  bo  positive  to  be  physically  meaningful  (and 

TE  tE 

as  we  shall  see  T.,  - A^.^o  and  T„  - A„  > 0),  then  we  require  that 

N N N N 


> 


for  a leading  edge  hinge 


or 


Oo  < 


for  a trailing  edge  hinge 
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in  order  to  be  able  to  achieve  a condition  of  zero  wall  interference 

7»  1 

by  increasing  the  non-dimensional  dynamic  pressure  from  0 to  A ' , 
MULTIPLE  WALL  SEGMENTS: 

Here  wo  briefly  indicate  the  generalization  to  multiple 
wall  segments  which  would  be  required  for  an  airfoil  with 
arbitrary  downwash.  The  basic  concept  is  that  each  segment 
has  a certain  angle  which  is  maintained  relative  to  those 
of  other  segments  by  feedback  control.  In  the  limit  of  a 
large  number  of  segments,  one  may  think  t>f  the  wall  having 
a continuous  variation.  Hence  the  generalization  of  (13) 
is 


w = - U a f ( X ) 

(23) 

-1^  - - 0 £ (X  - ‘'5-) 

where  f is  the  shdi e of  the  airfoil  (and  wind  tunnel  wall). 
Using  (23)  in  (10) , 


z=b/. 


P U 

' a>  03 


[a 


n=  1 


fix  + ^ 


- n(lb) 


- 0 f(x-^5)  ■ 20  ): 

n=l 


f (X 


- 


nhb)  ] 


(24) 


Using  (24)  in  (15)  again  leads  to  an  equation  of  the  form 
given  by  (16)  with  generalized  definitions  of  and  T^^. 
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RESULTS 


AERODYNAMIC  COEFFICIENTS: 

TE  TE 

T^  are  solely  functions  of  c/3b.  They 
are  plotted  in  Fig. 3.  For  convenience  they  are  re-normalized 
by  defining 


^N/(c/3b)^ 


(25) 


etc. 

For  large  etc.,  are  proportional  to  Hence 

for  large  c/^j^  or  small  3b/^,  approaches  a constant, 

namely  .5.  In  Fig. 4,  the  asymptotic  behavior  of  the  aero- 
dynamic coefficients  is  shown. 

AEROELASTIC  COEFFICIENTS  FOR  ZERO  WALL  INTERFERENCE: 

Pb 


The  variation  = ^^b^ ^ with 


for  leading  and 


trailing  edge  hinges  is  shown  in  Fig.  5.  Also  shown  for  re- 
‘V  C 2 

ference  is  = (^)  for  a trailing  edge  hinge.  Note 


T % 

that  to  insure  that  A < A 


pb 


. ''^TR 

, one  requires  6 , > A , 

'a  ' ;i<TE 

^N 


(-♦•  . 5 for  ~ 0)  . 


WALL  ROTATION  VS  AEROELASTIC  COEFFICIENT; 

In  Fig. 6,  the  wall  rotation  is  shown  as  a function  of 

various  0^^  . The  asymptotic  value  - *5 


N 


is  used.  Recall  equations  (18)  and  (19) . Note  that  0/^  = 0^^ 
at  “J:  = 0. 
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DISCUSSION  AND  CONCLUSIONS 

First  consider  the  rationale  for  a compliant  wall  wind 
tunnel . 

(1)  Any  movable  wind  tunnel  wall  will  be  inherently  com- 
pliant anyway,  so  that  one  must  take  this  into  account 
in  the  design  of  a nominally  rigid  wall. 

(2)  If  one  were  to  try  to  maintain  0/^  = 1 for  zero  wall 

interference  by  using  a rigid  wall,  one  would  have  to 

measure  6 and  a to  insure  that  6/  =1.  However,  for 

a 

a compliant  wall,  one  may  set  0 , appropriately  at 

''a 

A = 0 and  then  by  controlling  and  measuring  A insure 
that  0/^  = 1.  This  presumes  the  present  aeroelastic 
analysis  is  sufficiently  accurate,  of  course. 

The  principal  conclusions  from  the  present  analysis 

are 

(3)  The  aerodynamic  coefficients  associated  with  the 

compliant  wall  loading  rapidly  approach  their  asymp- 

fib 

totic  values  for  — < .5. 

c 

(4)  The  wall  segment  with  a leading  edge  hinge  is  more 
promising  than  one  with  a trailing  edge  hinge.  The 
latter  is  subject  to  aeroelastic  divergence.  Also  the 
compliant  wall  rotation  varies  more  rapidly  with  the 
aeroelastic  coefficient  about  its  optimum  value  for 
zero  wall  interference  for  a trailing  edge  hinge. 

Overall  the  compliant  wall  wind  tunnel  shows  real  promise 
as  a concept  for  reducing  wind  tunnel  wall  interference. 


even 
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In  closing,  it  should  be  noted,  that,  in  principle,  a i 

similar  concept  using  a porous  wall  may  be  employed.  This  con-  . i 

cept  is  briefly  explored  in  the  Appendix.  The  difficulty  with  | 

a porous  wall  is  the  basic  uncertainty  regarding  its  fluid 
mechanical  beljavior. 

i 


i 


I 
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AIRFOIL  IN  A WIND  TUNNEL 
FIG.  I 


SPRING  WITH 


COMPLIANT  WALL  GEOMETRY 


FIG.  2 


AERODYNAMIC  COEFFICIENTS  vs.  AIRFOIL  CHORD  TO  WIND  TUNNEL  HEIGHT  RATIO 


Pb/c 

ASYMPTOTIC  BEHAVIOR  OF  AERODYNAMIC  COEFFICIENTS  FOR 
SMALL  WIND  TUNNEL  TO  AIRFOIL  CHORD  RATIO 


AEROELASTIC  COEFFICIENTS  REQUIRED  FOR  ZERO  WALL  INTERFERENCE; 
DIVERGENCE  AEROELASTIC  COEFFICIENT  SHOWN  FOR  REFERENCE 


r 
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APPENDIX 


Porous  Wall  Wind  Tunnel 
The  field  equation  is 


3 d)  - (i)  = 0 

^ ^xx  ^zz 


The  boundary  condition  on  the  airfoil  is 


(A-1) 


z=0 


= w 


(A-2) 


where  w is  the  airfoil  downwash.  The  traditional 
boundary  condition  on  the  porous  wind  tunnel  wall  is 


P z 


z=b/. 


+ p u 4> 
* 00  00 


= 0 


(A-3) 


z=b/- 


where  b is  the  wind  tunnel  height  and  K characterizes  po- 

P 

rosity  of  the  wind  tunnel  wall. 

Again  using  a Laplace  Transform  with  respect  to  x, 
solving  the  resulting  ordinanry  differential  equation  in 
z and  using  the  transformed  boundary  conditions  gives  for 
the  Laplace  transform  of  the  velocity  potential,  (|)*, 


<P’ 


= - »*  in-f  I T"e-''®P'= 

z=0  n=0 


(A-4) 


where 


Y 5 [3  - P°o^°o]/[3+P°°^co] 


K 


K 


(A-5) 

P P 

Our  primary  interest  is  in  the  fluid  pressure,  p,  which  is 
determined  from  Bernoulli's  equation 


p = - p u (b 

^ f^oo  00 


(A- 6) 
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Using  (A-4)  , (A-6)  and  ‘Standard  inversion  formulae,  one 

obtains 


P U 

ss  ~ cm  r 


L -y  w(x-nBb) 
n=0 


+ Y.  Y (n+1)  pb) 


josing  y = 0,  (A-7)  reduces  to 


P U 

P = w(x) 


(A-7) 


(A-8) 


the  well  known  Ackeret  result  corresponding  to  zero  wall 
interference. 

From  (A- 5) , Y = 0 implies 


P U 

* <30  OO 


(A-9) 


Thus  by  varying  the  wind  tunnel  density,  for  example,  one 
could  ensure  that  (A-9)  is  satisfied  and  zero  wall  inter- 
ference is  achieved.  The  crucial  question  is,  to  what 
degree  does  (A-3)  accurately  describe  the  fluid  mechanics 
of  a porous  wall?  It  is  this  uncertainty  regarding  the 
adequacy  of  (A-3)  which  makes  the  compliant  wind  tunnel 
wall  concept  an  attractive  alternative. 

The  above  result  for  a porous  wall  wind  tunnel  has 
been  discussed  by  Maeder  and  Wood  and  was  well  known  to 
earlier  workers  on  wind  tunnel  wall  interference. 
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A Model  for  tho  Flow  tlu’oii[;h  Slots  tn  Wind  Tunnol  Walls 
llonuld  B.  Bliss,  I’rlnootou  Uoivorsity 

Introduot ion 

Ventilated  walls  are  useil  in  the  test  sections  of  transonic  wind 
tunnels  in  order  to  reiliice  hlocka^e  effects.  In  some  desi>;ns,  flow  is  allowed 
to  leave  or  enter  the  test  section  through  slender  slots  (larallel  to  the  flow 
wl\ich  are  distributed  alonj’.  the  test  section  walls.  I'he  si:e,  numlier,  aspect 
ratio  and  otl>er  geometric  properties  of  the  slots  may  vary  coi\siderahly  for 
different  tunnels.  In  order  to  include  wind  tunnel  wall  interference  effects 
in  transonic  flow  calculations  the  effect  of  the  sh>ts  on  the  wall  boundary 
condition  must  be  included.  A first  stej’  in  this  direction  is  to  understand 
the  flow  field  in  the  vicinity  of  a single  slot. 

The  type  of  slot  to  be  cotisidered  is  illustratetl  in  figure  I,  wltere 
the  flow  is  shown  leaving  the  test  section.  An  experimental  |>lot  of  pressure 
difference  across  a slot  in  uniform  flow  as  a function  of  the  mean  velocity 
through  the  slot,  suitably  nond imens iona 1 i ;ed , is  shown  in  I'ig.  J idoethert, 
1907).  Wheit  the  mean  velocity  is  small  the  behavior  is  linear,  whereas 
when  me  mean  velocity  is  large  the  behavior  is  more  nearly  ipiailratic.  If 
the  slot  were  subjected  to  a ui'imniform  flow  over  its  lengtlt  tlte  beliavivu' 
of  the  curve  would  be  altered  in  a way  related  to  the  streamline  curvature. 

The  fact  that  the  data  for  different  Mach  numbers  collapses  well  onto  a 
single  curve  suggests  that  the  behavior  of  the  slot  is  dominated  by  the  sub- 
sonic cross-flow  (y-z  plane)  and  indicates  that  the  slot  can  be  analyzed 
using  s lender-body  theory. 

Niunerous  investigators  have  studied  the  flow  through  slots 
annlytically  and  experimentally  ((loethert  , 19M;  Berndt  197.S,  1977).  There 


is  general  agreement  that  the  quadratic  behavior  seen  at  higher  values  of 
mean  slot  velocity  in  Fig.  2 is  associated  with  the  head  loss  as  the  cross- 
flow  separates  at  some  point  in  its  passage  through  the  slot.  The  dependence 
of  pressure  drop  on  the  square  of  velocity  is  the  same  as  that  encountered 
for  a jet  or  an  orifice  flow.  Perhaps  surprisingly,  the  linear  behavior 
that  occurs  when  the  mean  slot  velocity  is  small  has  not  been  properly 
explained.  Attempts  have  been  made  to  ascribe  the  linear  region  to  a viscous 
effect  in  the  vicinity  of  the  slot.  This  explanation  appears  to  be  incorrect; 
in  fact  the  effect  of  viscosity  is  probably  to  contribute  another  quadratic 
term.  Another  possibility  is  that  the  linear  term  is  related  to  an  inter- 
action of  the  flow  with  the  trailing  edge  of  the  slot.  For  instance,  a 
blunt  trailing  edge  would  cause  a flow  deflection  not  unlike  that  produced 
by  the  holes  in  a perforated  wall  which  are  known  to  have  a linear  character- 
istic. However,  if  this  were  the  mechanism,  the  slope  of  the  linear  region 
would  probably  be  dependent  on  Mach  number.  Also,  both  sharp  and  blunt 
edged  slots  are  known  to  exhibit  linear  behavior  at  low  mean  slot  velocities. 

In  the  analysis  that  follows,  the  linear  behavior  is  shown  to  be 
a consequence  of  the  correct  application  of  slender  body  theory  to  this 
problem.  Physically,  the  linear  behavior  is  associated  with  the  region 
downstream  of  the  leading  edge  where  the  free  surface  flow  is  beginning  to 
form  but  has  not  reached  its  final  cross-flow  configuration  as  an  orifice  or 
jet  flow.  If  the  mean  flow  through  the  slot  is  small,  the  cross-flow  never 
reaches  its  fully  developed  state  and  the  linear  term  dominates.  When  the 
mean  flow  is  large  most  of  the  cross-flow  resembles  a fully  developed  jet  or 
orifice  flow  and  the  quadratic  term  dominates.  The  somewhat  analogous  situation 
in  classical  slender  body  theory  is  that  lift  forces  arise  only  from  regions 


in  which  the  cross-flow  experiences  a streamwise  rate  of  chanRc,  e.g.  only 
from  regions  of  the  body  where  the  cross-sectional  area  is  changing. 

In  addition  to  showing  how  the  linear  behavior  arises,  the  analysis 
also  shows  how  to  collapse  the  data  in  order  to  compare  the  features  of 
slots  of  different  aspect  ratios. 


Formulation 


For  reasons  of  simplicity  oitly,  the  discussion  is  restricted  to 
incompressible  flow.  Assume  that  the  fluid  is  inviscid  and  irrotational , 
then  the  flow  is  governed  by  Laplact''s  equation 


0 


where  ? (x,y,2)  is  the  velocity  potential  such  that  the  velocity  vector  is 


given  by 


The  configuration  of  the  slot  and  the  coordinate  system  are  as  shown  in  Fig. 
1.  The  slot  is  assumed  to  lie  in  an  infinite  plane  below  a uniform  free 
stream. 

The  problem  will  be  formulated  in  terms  of  the  method  of  matclied 
asymptotic  expansions.  The  velocity  potential  can  be  expressed  as 

Where  ^ is  a nondimen sional  perturbation  velocity  potential.  Notice  that 
when  there  is  no  pressure  difference  across  the  plane  containing  the  slot, 
i.e.  Ap  — 0 , then  the  effect  of  the  slot  vanishes.  The  inner  length 


scale  is  defined  to  be  the  slot  width  a.  The  relevant  small  parameter  for 


which  is  a measure  of  the  flow  deflection  angle  in  the  slot.  Then  the  outer 
length  scale  is  defined  to  be  L - Ye  . Notice  that  the  outer  length 
scale  is  not  the  slot  length,  , which  is  not  really  the  appropriate  scale 
for  variations  in  the  streamwise  direction.  In  fact,  this  procedure  could 
be  used  to  solve  the  problem  of  a semi-infinite  slot  in  which  there  is  no 
possibility  of  using  slot  length  as  an  outer  length  scale. 

The  following  nondimensionalization  is  used  for  the  variables: 

1 = Llfci* 

u= 

K=  U*,r 

This  leads  to  a non-dimensional  outer  expansion  of  the  form: 

i= 

Substituting  into  Laplace's  equation  gives 

The  inner  variables  are  defined  by  a stretching  of  the  cross-flow 
coordinates 


(note  that  y = y/a,  2 = z/a) 

Then  assume  an  inner  expansion  of  the  form 


Substituting  into  Laplaces  equation  gives,  to  lowest  order 

which  shows  that  the  inner  region  behaves  as  a two  dimensional  potential  flow 
in  the  cross  flow  plane.  Solutions  must  be  of  the  form 

The  boundary  conditions  for  the  problem  are  that 


OS,  ^/=r  X V 


^0 

7=0 


on  tlae  ricjict  wall 


p - U'Asi&vxf  r ovi  ^urfac^  i.t  Ov  ilof 


The  first  and  second  of  these  conditions  are  applied  to  the  outer  velocity 
potential  fkf.zi  , whereas  the  second  and  third  conditions  are  applied 
to  the  inner  velocity  potential  . Any  indeterminancy  that 

remains  after  the  application  of  the  appropriate  boundary  conditions  to  each 
solution  is  resolved  when  the  solutions  are  matched.  For  the  outer  solution, 
to  the  present  order  of  expansion,  the  slot  appears  as  a line  of  sinks  dis- 
tributed along  the  x-axis.  The  strength  of  this  distribution  is  determined 
by  matching  with  the  inner  solution,  as  is  usual  with  slender  body  theory. 

The  boundary  conditions  for  the  inner  solution  require  closer  examination. 

Bernoulli's  equation  in  dimensional  form  is 


After  nondimensionalizing,  re-expressing  in  inner  variables,  and  setting 
p * Pj  this  equation  becomes 


Retaining  only  the  lowest  order  terms  gives 


L St" 




This  is  the  condition  that  the  free  surface  be  at  constant  pressure.  The 
unknown  position  of  tlie  free  surface  expressed  in  inner  variables  is 


The  solution  must  also  satisfy  the  condition  that  the  flow  be  tangent  to 
the  free-surface  at  the  free-surface  location.  If  the  free  surface  is  des- 


cribed by 


then  in  dimensional  coordinates  the  flow  tangency  condition  is 

=0 

Re-expressed  in  inner  variables  this  becomes 

aS  + i.  al  5i  . I 3B  aj  _ n 

a?  5| 

Letting  ^ using  the  expansion  for  5.  gives 

afif  + i.?f  =0 


Rctflining  the  lowest  order  terms  gives  the  flow  tungeiicy  condition  in  inner 


variables 


-^21^ 

1.  an*  ay  du 
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Solutions  and  Matching 


The  general  free  surface  problem  in  the  slot  would  be  extremely 
difficult  to  solve.  Therefore  this  section  is  restricted  to  general  comments 
on  the  solution  structure  and  on  the  matching  procedure.  As  stated  before 
the  inner  solution  must  have  the  form 


As  is  usual  for  slender  body  theory,  far  from  the  slot  the  velocity  potential 
will  look  like  a simple  source  or  sink  whose  strength  is  related 
to  the  rate  of  change  of  cross-sectional  area  of  the  flow  in  the  slot. 

+q(r}sf° 


Where  r - j and  is  the  outer  limit  of  the  inner  solution.  Thus 

.where 

equals  the  cross-sectional  area  occupied  by  the  flow  in  the  slot  (dimen- 
sionally S - Q*  S ) . 


The  outer  solution  can  be  expressed  as  an  inteRral  of  a source- 


sink  distribution  over  the  slot 


+ -L- 

^ zif  JjiAv* 


where  1 - t/i^  is  the  nondimcnsional  slot  length. 

This  integral  can  be  shown  to  have  the  following  fonn  when 

is  small  (see  Ashley  and  Landahl,  1965) 


Re-expressing  the  outer  solution  in  inner  variables  and  using  the  expression 
for  small  gives  the  inner  limit  of  the  outer  solution. 

o,*‘-  jv'uF  _ + X r'r(c)-nx‘) 

'f  =-  ^ 2Tr  Tx--ci  ''  ' 

Here  /olG  is  treated  as  an  Oil]  quantity,  as  is  sometimes  done  in  slender- 
body  theory.  Alternatively,  this  problem  can  be  treated  more  formally  by  a 
modification  of  the  expansion  procedure,  but  the  results  would  be  the  same. 

The  matching  procedure  requires  that  the  inner  limit  of  the  outer 
solution  equal  the  outer  limit  of  the  inner  solution:^  -y  . The  match- 


ing shows  that 


S(<*)  = ?(<♦) 


and  then 


_ -Sid . _eL__  + X r %)_- ikl)  j y- 

Now  the  problem  for  flow  in  the  inner  region  can  be  restated.  The 
velocity  potential  is  of  the  form 


q>  = j + ft*) 

dij  01" 

Far  from  the  slot  Cn  t~  ^ -^OO, 

In  the  slot  the  constant  pressure  boundary  condition  becomes 


' \Ooj 


\ 3=  l-^</UV 


The  flow  tangency  condition  is 


^i}rj  ^ 

These  equations  constitute  the  final  form  of  the  inner  problem.  In  addition, 
depending  on  the  assumed  shape  of  the  free  surface,  other  conditions  may  be 
required  at  the  slot  edges.  For  instance,  the  sides  of  a sharp  edged  slot 
may  require  a Kutta  condition.  The  difficulty  in  solving  the  problem  as 
posed  above  is  that  both  the  free  surface  shape  and  the  velocity  potential 
are  unknown.  Certain  limiting  cases  are  solved  below. 


Slot  Flow  for  Small  Free  Surface  Displacements 


V*<  4 4 

^ ^ . Under  these  conditions,  since  € measures  the  flow  deflec- 


,♦ f ^ 

The  nondimensional  slot  length  iS  1 . When  t # then 

A 


tion  angle  in  the  slot,  the  free  surface  displacement  is  sufficiently  small 


the  quadratic  terms  in  the  pressure  boundary  condition  can  be  nejilected 
(they  are  of  order  ^ );thus  ^ • 

The  followinyi  function  will  satisfy  this  simplified  pressure  boundary  condition; 


where  w is  a constant  of  order  unity  (to  be  discussed  later).  Differentiating 
and  substituting  for  g(x*)  yields 

f* 

- S(0-x  H-C,  + ^4;^  -zf  ''  ' 


This  integral  equation  must  be  solved  for  S'(x*).  The  presence  of  the 
arbitrary  constant  should  allow  the  application  of  a Kutta  condition  at 
the  front  of  the  slot,  namely  S' (0)  = 0* 

Now  that  the  equation  has  been  formulated,  it  is  convenient  to  change  to 
a new  dimensionless  variable  Si  = x/!l,  so  that 


* f ^ 

X = ej-x 

Cl 


$'(<*)=  S'(x)|.  = s\i)4i 


and  the  integral  equation  becomes 


where  Cj^  is  a redefinition  of  the  arbitrary  constant. 

The  velocity  potential  which  satisfies  the  boundary  conditions  on  the 
crossflow  is  given  by 

{+  M-2?-^ATu!rl  - 1] 


I 


is  a complex  quantity.  This  velocity  potential  is  asso- 


where  ^ — 2 -V 

ciated  with  the  potential  flow  through  a slit;  the  velocity  field  exhibits 
a square-root  singularity  at  the  side  edges.  A brief  calculation  gives 


The  integral  equation  now  becomes 


Clearly,  from  the  form  of  this  integral  equation,  the  solution  must 
have  the  form 

Recall  the  dimensional  relationships;  %-\  X=1  fihici  ) = SLx) . 

Thus 

where  and  the  function  G is  defined  by  the  above  expression. 

T 

The  volume  flow  rate  from  the  slot  is,  to  lowest  order, 

Q = U„;S?  ^ur^a  i 

where  is  the  mean  velocity  through  the  slot.  Substituting  the  above 
expression  for  S gives 


or 


This  shows  a 1 inoar  relationship  between  the  pressure  difference  across  the 
slot  and  the  mean  velocity  through  the  slot.  Recall  that  this  rclationshi[i 
will  hold  only  for  small  free  surface  displacements. 

The  relationship  between  pressure  difference  and  mean  slot  velocity  can  also 
be  written  as 

(f  a - 

Again,  this  linear  relationship  holds  only  for  » which  corres- 

' l/oo 

ponds  physically  to  small  free  surface  di siilacements  compared  to  the  slot 
width. 

The  slope  of  the  pressure  difference  versus  velocity  curve  is 

The  quantity  a/£  is  the  aspect  ratio  of  the  slot.  Its  appearance  is  analogous 
to  the  familiar  result  in  slender  wing  theory  that  the  lift  curve  slope  is 
proportional  to  the  aspect  ratio.  The  fact  that  the  argument  of  the  function 
G is  logarthmic  suggests  that  its  dependence  on  slot  aspect  ratio  may  he 
relatively  weak.  The  integral  equation  must  bo  solved  to  determine  the  func- 
tion 

In  order  to  estimate  the  type  of  results  given  by  this  analysis,  an 
approximate  solution  to  the  integral  equation  is  now  obtained.  For  very  small 
aspect  ratio,  , the  term  involving  sX^)  which  has  the  coefficient 
will  dominate.  Then  the  equation  is  approximately 


j 


rhysically,  this  ov\uation  halaucos  tl\o  aocolorat  ion  of  fluid  noai*  tho  slot  ami 
the  ai'plied  pressvwe  d i f ferenoe  Ap  ; tl\o  effect  of  the  source-sink  distrlhu- 
tion  on  the  juessure  distrihution  in  the  slot  is  lu'nlected.  iActvially  tliis 
approximation  is  not  formally  consistent  since  it  was  previously  assumed  that 
e«,  <v_(  and  ivit  = cyii|  in  the  derivation  of  the  etpiation.l  Setting; 

C,-0  , so  that  S (O)  - 0 . ^ives 


Inte^rat  inj;  >;ives 


from  which 


4 

cv 


This  result  >;ives  a fuitction  lij  which  can  l>e  viewed  as  a first  approximation 


to  the  actual  funct  ion  (I: 


As  a refinement,  the  nej;lected  terms  can  he  include»l  hy  approximating 

as  a linear  function,  a form  sug^’ested  hy  the  approximate  solution  above. 
One  possibility  is  simply  to  use  this  approximate  solution  in  the  terms 
previously  neijlected.  Taking  this  apj'roach  and  suhst  i t ut  inn  into  the  inteijral 
equation  Rives 


After  choosing 
becomes 


to  satisfy  the  Kutta  condition,  the  equation 


I 


(K 


Although  the  expression  for  Sit)  is  now  singular  at  l(  = ( , the  function 

is  integrablo.  The  result  is 


Thus  the  function  G2  which  is  the  second  approximation  for  G is 


81 


For  a slot  of  aspect  ratio  - 0.  i , the  resulting  values  of  G^  and 


G2  are 


(Vj  S'.iff  2S  4. 

The  corresponding  relationship  between  pressure  difference  and  mass  flow  rate 
are 


I 


0»S(o 


These  numbers  are  considered  to  be  in  qualitative  agreement  with  the  avail- 
able experimental  data  (e.g.  see  Fig.  2). where  the  slope  of  the  pressure  versiJs 
flow  rate  curve  seems  to  fall  in  the  range  of  less  than  0.5  to  greater  than 
unity  (Goethort,  1961).  Aside  from  the  approximate  nature  of  the  solution 
presented  here,  there  are  several  other  factors  which  could  affect  the  accuracy 
of  the  calculation.  The  presence  of  a wall  boundary  layer  and  a shear  layer 
in  the  slot  could  substantially  lower  the  magnitude  of  the  velocity  in  the  slot. 


Since  transonic  tunnel  speeds  are  of  primary  interest,  the  effect  of  com- 
pressibility must  also  be  considered.  All  of  those  effects  are  believed  to 
act  to  increase  the  slope  of  the  pressure  drop  versus  flow  rate  curve.  Regard- 
less of  these  refinements,  the  analysis  demonstrates  that  there  is  a (pre- 
viously unrecognized)  mechanism  for  linear  behavior  of  this  curve  within  the 
framework  of  an  inviscid  slender  body  theory  analysis. 

Fully  Developed  Flow  Through  the  Slot 


When  1 the  cross-flow  in  the  slot  will  closely  resemble 

a fully  developed  jet  or  orifice  flow  (depending  on  the  slot  cross-sectional 

geometry).  In  this  section,  the  problem  fonnulation  presented  previously  is 

shown  to  be  consistent  with  this  picture  of  the  flow  field.  In  this  case 

1 • Away  from  the  ends  of  the  slot  the  velocity  potential 

^ a ^ 

is  expected  to  be  independent  of  fo  lowest  order  and  ' 

the  constant  local  flow  rate  through  the  slot.  Under  these  conditions. 


F' 
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Completely  neglecting  the  integral  that  appears  in  the  exact  expression  for 
means  that  end  effects  are  being  ignored.  Differentiating  the 


above  gives 


This  quantity  is 


01 W 


away  from  the  ends  of  the  slot  and  can  be  neglected 


d<l> 

in  the  constant  pressure  boundary  condition.  The  quantity  can  also  be 

neglected  by  the  (consistent)  assumption  that  the  velocity  potential  is 
independent  of  to  lowest  order.  Therefore,  the  pressure  boundary  condition 


becomes 


UlfJ 


Thus,  to  good  approximation,  the  velocity  potential  for  the  equivalent 
steady  two-dimensional  free  surface  flow  can  be  used,  provided  the  changes 
in  free  surface  shape  are  confined  to  the  region  where  the  flow  has  become 
essentially  parallel  and  the  pressure  is  nearly  constant  (=  p^).  This 
situation  is  illustrated  schematically  in  Figure  3 for  the  case  of  a sharp 
edged  slot. 

To  obtain  the  relationship  between  pressure  difference  and  flow 
through  the  slot  it  is  not  necessary  to  consider  the  velocity  potential  in 
any  detail.  Suppose  that  the  final  width  of  the  jet  leaving  a slot  of  width 
unity  is  CT'  . The  precise  value  of  o'*  depends  on  the  slot  geometry,  but 
typically  o'  - 0.611  for  a sharp  edged  slot,  and  0"^=  1.0  for  a thick  slot 
without  internal  separation.  Applying  the  boundary  condition  to  the  cross 
flow  in  the  region  downstream  of  the  slot  where  the  flow  is  almost  parallel 


gives 


=1 


where  is  the  final  cross  flow  velocity  downstream  of  the  slot.  Then 


Integrating, 


S'(x')  = S(  = = s' 

S(r)=(j'x* 


A constant  of  integration  that  could  be  included  in  the  above  has  been 
omitted  since  end  effects  are  being  neglected.  Set  x"*-  ('*  and  re-express  in 
dimensional  form  t 


— - ^ 


As  before,  the  flow  rate  through  the  slot  is 


(?  = (XS  = ( 


Hence 


Voo{^aicr^  = ur^(Ki 


or 


Ap L- 

■7*  — at'Z- 
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For  a sharp  edge  orifice  ~ 2.68  • This  equation  was  used  to  generate 

the  "quadratic  behavior"  curve  in  Fig.  2 assuming  a sharp  edged  orifice;  note 
that  the  agreement  with  the  data  is  quite  good. 

Note  that  the  above  result  can  also  be  rewritten  as 
.z\  / . > . \2- 
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Since  the  equation  for  linear  behavior  could  also  be  written  in  an  analagous 
form,  it  is  suggested  that  all  slender  slot  behavior  can  probably  be  expressed 
as 


In  fact  this  can  be  seen  directly  from  the  final  statement  of  the  inner  pro- 
blem. For  the  above  form  to  hold,  the  cross-sectional  shape  of  the  slots 
must  be  geometrically  similar  and  there  must  be  no  boundary  layer  effects. 
The  parameter  UL  may  prove  to  have  only  a weak  influence.  It  would  be 
most  interesting  to  plot  an  appropriate  set  of  experimental  data  in  the  form 

versus 

and  see  to  what  extent  it  collapses  onto  a single  curve. 
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